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We study finite size effects in the variance of the displacement of a tagged particle in the station- 
ary state of the Asymmetric Simple Exclusion Process (ASEP) on a ring of size L. The process 
involves hard core particles undergoing stochastic driven dynamics on a lattice. The variance of 
the displacement of the tagged particle, averaged with respect to an initial stationary ensemble and 
stochastic evolution, grows linearly with time at both small and very large times. We find that at 
intermediate times, it shows oscillations with a well defined size-dependent period. These oscilla- 
tions arise from sliding density fluctuations (SDF) in the stationary state with respect to the drift 
of the tagged particle, the density fluctuations being transported through the system by kinematic 
waves. In the general context of driven diffusive systems, both the Edwards- Wilkinson (EW) and 
the Kardar-Parisi-Zhang (KPZ) fixed points are unstable with respect to the SDF fixed point, a 
flow towards which is generated on adding a gradient term to the EW and the KPZ time-evolution 
equation. We also study tagged particle correlations for a fixed initial configuration, drawn from 
the stationary ensemble, following earlier work by van Beijeren. We find that the time dependence 
of this correlation is determined by the dissipation of the density fluctuations. We show that an 
exactly solvable linearized model captures the essential qualitative features seen in the finite size 
effects of the tagged particle correlations in the ASEP. Moreover, this linearized model also provides 
an exact coarse-grained description of two other microscopic models. 

PACS numbers: 02.50.Ey, 05.40.-a, 05.60.Cd, 66.30.Dn, 05.50.+q 



I. INTRODUCTION 

The Asymmetric Simple Exclusion Process (ASEP) is 
a prototype model to study driven diffusive motion 
The process involves hard-core particles hopping between 
neighboring sites of a lattice, with an asymmetry in hop- 
ping rates which incorporates the effect of an external 
drive. In the case of completely asymmetric dynam- 
ics, the ASEP reduces to the totally asymmetric simple 
exclusion process (TASEP). In either case, the system, 
on a periodic lattice, settles at long times into a non- 
equilibrium stationary state in which each particle has 
an average drift velocity, and there is a finite current 
through the system. 

In this paper, we are interested in the temporal growth 
of the fluctuations in the displacement of a tagged ASEP 
particle in the steady state in a finite system of size L 
in one dimension. We monitor the variance a 2 (L,t) of 
a tagged particle around its average in time t, starting 
from the stationary ensemble of configurations. Some 
asymptotic results are already known 0, d, @] ■ These 
are summarized below. 

(a) On taking the limit L — > oo first, followed by the 
limit t — > oo, the fluctuations are diffusive, growing lin- 
early in time: a 2 (L,t) ~ D^t. Here, Do is a known 
function of the particle density and the external bias [2| , 
and is given later in the paper. 
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FIG. 1: (Color online) Monte Carlo (MC) simulation re- 
sults for the time dependence of the variance a 2 {L, i) for the 
TASEP on a finite lattice, averaged over 10 5 MC runs. The 
different time regimes are also marked. Here, L = 1024, the 
particle density p — 0.25. 



(b) In the opposite limit with t — > oo before L — > oo, 
a 2 (L,t) again behaves diffusively, i.e., cr 2 (L,t) ~ D(L)t, 
but now the diffusion constant D(L) depends on the sys- 



2 



tern size L and scales as D(L) ~ for large L Q. 
A natural question arises: What is the full behavior of 
a 2 (L,t) as a function of time in between these two ex- 
tremes? In this paper, we show that in the intermediate 
regime, <r 2 (L,t) shows striking oscillations as a function 
of time. Indeed, there are two time-scales in the problem, 
Ti ~ L and T 2 ~ L 3 / 2 . The limiting behaviors (a) and 
(b) hold in the regimes t <C T\ and t ^> T2, respectively, 
and in between (T\ <C i -C T 2 ), the quantity a 2 (L,t) 
oscillates with time, as shown in Fig. [TJ 

In measuring the fluctuations, if one does not start 
from the stationary ensemble, but instead from an arbi- 
trary but fixed configuration drawn from the stationary 
ensemble, the corresponding variance s 2 (L,t) of the dis- 
placement behaves very differently from a 2 (L,t). The 
behavior of s 2 (L,t), in the limit of an infinite system, is 



known: s (t) = lim^_ 



! (T,i) ~ t 2 ' 3 



For a finite 



system, we find a similar growth until a characteristic 
time T* ~ L 3 / 2 , after which this correlation grows lin- 
early in time: s 2 (L,t) ~ D{L)t with D(L) ~ 

Some of the results for the tagged particle correlations 
discussed in this paper were known previously. However, 
a complete and a consistent picture on a finite lattice is 
lacking in the literature. The principal purpose of this 
paper is to discuss finite size effects on the tagged particle 
correlation in the ASEP on a one-dimensional ring from 
a unified point of view. In order to do this, we include 
some known results for completeness, in addition to the 
new results on the finite size effects. The main results 
are summarized below. 

1. The behavior of er 2 (T, t) as a function of time can be 
characterized by two time scales T\ ~ L and T 2 ~ T 3 / 2 , 
as shown in Fig. [1] 

(a) In the initial linear regime (t <C Ti), the variance 
a 2 (L,t)~D t 0. 

(b) In the oscillatory regime (Ti <C t <C T2), the quan- 
tity cr 2 (L,t) oscillates as a function of t. The amplitude 
of oscillations is proportional to the system size L while 
the time period of oscillations is T = L/u, where u is the 
mean velocity of the stationary state density fluctuations 
relative to the average drift velocity of the particles. 

(c) In the late time regime (t 3> T2), the quantity 
a 2 (L,t) ~ D{L)t with TJ>(L) - 

2. The time variation of s 2 (L,t) can be characterized 
by a single time scale T* ~ L 3 ! 2 . The variation of s 2 (L, t) 
with time for two different system sizes is shown in Fig. 
® 

(a) t < T*: In this regime, s 2 {L, t) ~ t 2/3 . 

(b) t > T*. Here, s 2 (T,i) ~ D(L)t with D(L) ~ -±=. 

The long-time behavior of both the functions a 2 (L,t) 
and s 2 (L, t) is diffusive, with the same diffusion constant 
D(L) which scales with the system size as D(L) ~ -j=. 
This behavior is attributed to the motion of the center- 
of-mass. 

The physical origin of oscillations in a 2 (L,t) is the 
occurrence of kinematic waves that transport stationary 
state density fluctuations through the system (Appendix 



[5]) . Kinematic waves are known to arise in a variety 
of circumstances involving flow e.g., in flood movement 
in long rivers ||, traffic [2, [nj, flow of granular parti- 
cles through vertical tubes and ho ppe rs 111 ] . motion of 
transverse fluctuations in interfaces [l2,[l3|, field- induced 
transport in random media as in drop-push dynamics 
[Til ]. In our system, the velocity of the sliding density 
fluctuations (SDF) relative to the average drift of the 
ASEP particle determines the period of the oscillations. 
Eventually, the oscillations damp down as the density 
fluctuations dissipate owing to the stochasticity and non- 
linearity inherent in the dynamics. This dissipation in 
time is well captured by the temporal behavior of the 
function s 2 (T, t), as we show in this paper. 



3. The ASEP can be mapped to a non-equilibrium 
model of interface growth in the Kardar-Parisi-Zhang 
(KPZ) universality class [l5|. The resultant time evolu- 
tion equation for the interface is the usual KPZ equation, 
augmented by a drift term which accounts for the SDF. 
In this paper, we consider the interface equation with 
the drift term in the linear approximation, and solve ex- 
actly for a 2 (L,t) and s 2 (L,t). Our analytical solution 
of the linearized model captures the essential qualita- 
tive features seen in these two quantities for the ASEP. 
Moreover, we show that the linearized model is an ex- 
act coarse-grained description of two specific microscopic 
models of interacting particles, the Katz-Lebowitz-Spohn 
(KLS) model [l6| at a specific value of the temperature 
and the Asymmetric Random Average Process (ARAP) 
[l7| . Our analytic solution quantitatively describes the 
behavior of the tagged particle correlations for these two 
models. 



The paper is organized as follows. In Section [TTJ we 
define the model and the variances a 2 (L,t) and s 2 (L,t) 
of the tagged particle process, followed by a discussion in 
the context of the KPZ interface equivalent to the ASEP. 
In Section llTTI we give a representation for the tagged par- 
ticle displacement on an infinite lattice, and show how it 
correctly accounts for the known results for the variance. 
This is followed by Section IIVI where we adduce physi- 
cal arguments for the observed behavior of a 2 (L,t) and 
s 2 (L,t), both on an infinite as well as a finite lattice. 
The time-evolution equation for the interface-equivalent 
of the ASEP is solved exactly in the linear approximation 
in Section fVl which explains some qualitative features of 
the tagged particle correlations for the ASEP. As we show 
in Section [VI 1 the exact solution also explains quantita- 
tively the variance of the tagged particle displacement 
for the KLS model at a specific value of the temperature 
and for the ARAP. In Section IVlIl we discuss the center- 
of-mass motion for the ASEP, followed by conclusions in 
Section E!IE1 
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FIG. 2: (Color online) MC simulation results for the time 
dependence of the variance of the displacement for the TASEP 
on a finite lattice for a 2 (L,t) and s 2 (L,t). In both cases, the 
averaging is over 10 MC runs. Here, the density p = 0.25. 
The two system sizes are 512 and 1024. 



II. TAGGED PARTICLE CORRELATIONS 
A. The Model 

We consider the ASEP on a periodic lattice of L sites. 
N indistinguishable hard core particles are distributed 
over the lattice sites with each site either singly-occupied 
or empty. The particle density p = j- is held constant 
when the limit L — > oo, N — > oo is taken. Let rii = 0,1 
with i = 1,2, ... , L denote the occupancy of the i-th 
site. The system evolves according to a stochastic dy- 
namics: during an infinitesimal time interval dt, a parti- 
cle attempts to hop to the site to its right with probabil- 
ity pdt, to the left neighboring site with probability qdt, 
and continues to occupy the original site with probability 
1 — (p + q)dt. The attempted hop is successful only if the 
sought site is empty before the hop. L attempted hops 
constitute one time step. Clearly the total number of par- 



ticles is conserved under the dynamics. For the TASEP, 
the motion of the particles is entirely in one direction i.e., 
p = 1 and q = or, vice-versa. Note that with p = q, 
each particle moves symmetrically to the left and to the 
right. The model then reduces to the symmetric Simple 
Exclusion Process (SEP), an equilibrium model of hard 
core particles diffusing on a lattice. 

It is evident from the rule of evolution that the dy- 
namics is ergodic. In the limit t — > oo, the system settles 
into a nonequilibrium stationary state. Ergodicity en- 
sures uniqueness of the stationary state. On a lattice of 
size L, the stationary state is one in which all configu- 
rations C have the same weight p(C) = EILIlzIQI pj ft 
then follows that the two-point correlation function in 
the stationary state riiUj — xf^izjj for i =/= j- Here, as in 
the rest of the paper, the overbar is used to denote an av- 
erage over the stationary ensemble of configurations. In 
the thermodynamic limit, the steady state has a prod- 
uct measure form. In the stationary state, the ASEP 
supports a steady curren t of particles whose mean value 
J = (p — q)rii(\ — Uj+i) for a system of size L is given 
by J = (p — q)p(l — p) + 0(4). Correspondingly, the 
mean velocity of a particle in the stationary state, given 
by v P = j, equals v P = (p- q)(l - p) + 0(±). 

Besides particle motion, there is also a motion as- 
sociated with the coarse-grained density fluctuations in 
the stationary state of the ASEP. The density fluctua- 
tions are transported as a kinematic wave with velocity 
Vk = fp ( see Appendix fA"]) . Thus, relative to the aver- 
age drift of the particles, the density fluctuations 'slide' 
with velocity vk — Vp. We refer to this relative motion 
as the sliding density fluctuations (SDF). Stochasticity 
and nonlinearity in the dynamics lead to dissipation of 
the density profile so that the wave of fluctuations ulti- 
mately dies down in time. The kinematic wave velocity 
for an infinite system is given by vk = {p — <?)(1 — 2p). 
This velocity, derived on the basis of hydrodynamic con- 
siderations discussed in Appendix also appears as the 
imaginary part of the low- lying eigenvalues of the Markov 
matrix governing the time evolution of the ASEP in the 
totally asymmetric case q = [HI . Figure [3] shows the 
time evolution of a stationary ASEP configuration. One 
clearly sees that the space-time trajectory of a particle is 
distinctly different from that for a coarse-grained density 
fluctuation. 



B. Two different observables 

The primary quantity of interest in this paper is the 
variance of the displacement of a tagged particle around 
its mean trajectory. Note that the displacement y(n, t) 
of the n-th particle in time t ranges from — oo to oo even 
though the system size L is finite. 
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FIG. 3: (Color online) Time evolution for a stationary TASEP 
configuration. Here, the system size L = 200, the density 
p — 0.4. The inverse of the slope of the thick black line is 
the velocity of the tagged particle vp « (p — q)(l — p) = 0.6 
while that of the thin red line is the velocity of the density 
fluctuations vk ~ (p — g)(l — 2p) = 0.2. 

1. Average over histories and initial conditions: a 2 (L,t) 

Here we start with the stationary ensemble of the 
ASEP configurations at t = and monitor the motion of 
a tagged particle. Let A n (t) denote the deviation in the 
displacement of the n-th tagged particle from its average 
value in time t, the averaging being with respect to both 
the initial stationary ensemble at t = and stochasticity 
in the evolution of configurations. In symbols, 

A n (t) = y(n, t) - y(n, 0) - ([y(n, t) - y{n, 0)]) 

[Ensemble of Initial Conditions] . (1) 

The angular brackets denote averaging over the stochas- 
ticity in the evolution of configurations, while the over- 
bar is used to denote averaging with respect to the ini- 
tial stationary ensemble at t = 0. In the function A n (t), 
the quantity (y(n, t) — y(n, 0)) gives the average displace- 
ment of the n-th tagged particle in time t in the station- 
ary state, and hence is equal to vpt. The variance of the 
displacement, measuring the stationary state fluctuations 

I 



where the angular brackets denote averaging with respect 
to stochastic evolution. 

(1) On taking the limit to — ► 0, when the averaging 
is only with respect to stochastic evolution, we get the 
function s 2 (L, t). 



in A„(i), is given by 

a 2 (L,t)^jAjW)- (2) 

Since the system is translationally invariant, the above 
quantity is the same for all particles, and hence we do an 
additional averaging of a 2 (L,t) over the particles. 

In this paper, we study how the finiteness of the system 
size L affects the behavior of a 2 (L,t) in time. 



2. Average over histories only: s 2 (L,t) 

Beginning with an arbitrary but fixed configuration, 
chosen from the stationary ensemble and following its 
evolution in time, it is interesting to monitor the fluctua- 
tions in the displacement of a tagged particle, and study 
how finite size effects come into play here. These fluctu- 
ations were first studied by van Beijeren for an infinite 
system in Q. In this case, one averages only over the 
stochasticity in the evolution of the configuration. We 
define d n (t) to be the deviation in the displacement of 
the n-th tagged particle from its mean value in time t, 
starting from an arbitrary but fixed configuration, drawn 
from the stationary ensemble at t = 0. Explicitly, 

d n (t) - y(n, t) - y(n, 0) - {[y(n, t) - y{n, 0)]) 

[Fixed Initial Condition] . (3) 

The corresponding variance is given by 

s 2 (L,t)^(dl(t)). (4) 

The system being translationally invariant, the above 
quantity is also averaged over the particles. 

We note that the behavior of both a 2 (L, t) and s 2 (L, t) 
can be extracted from two different limiting values of a 
single function that measures the fluctuation in the dis- 
placement of a tagged particle in the stationary state. At 
time t = 0, we start with an arbitrary but fixed config- 
uration Co, drawn from the stationary ensemble of con- 
figurations. After Co has evolved in time for an interval 
to, we start measuring the variance Ct, (to, to + t) of the 
displacement of a tagged particle around its average: 



(5) 



(2) On the other hand, since the system is ergodic, in 
the limit to — > oo, the initial configuration Co evolves 
into an ensemble of stationary configurations. Thus, the 
function lim^^oo Cl(£o,£o + t) measures fluctuations in 
the stationary state where the averaging is with respect 



C L (to, h+t) = ([y(n, to + t)- y{n, t ) - ([y(n, t +t)- y{n, to)])] 2 ), 

I 
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to both the initial stationary ensemble and stochastic 
evolution. Hence, this quantity is identical to a 2 (L,t). 



C. Interface equivalent 

Many of our results for a 2 (L,t) and s 2 (L,t) can be 
explained in terms of the continuum description of the 
ASEP, achieved by mapping it to a model of nonequi- 
librium interface growth [5| that belongs to the KPZ 
universality class of noncquilibrium interface dynamics 
[la |. The mapping that we use in this paper employs 
the tagging process of the particles in a direct and es- 
sential way in the translation [lj|. Let the particles be 
labeled 1, 2, .., N sequentially at the initial instant. Since 
the particle motion is in one dimension, and there is no 
overtaking, this labeling will be preserved for all subse- 
quent times. The corresponding interface is obtained by 
identifying the tag label n with the horizontal coordinate 
n for the interface, while the set {y(n,t}}, denoting the 
particle locations at time t, maps onto the set of local 
interface heights {h(n,t)} defined by 



h(n,t) = y(n,t) - 



(6) 



Since the inequality y(n + > y(n,t) + 1 holds, it 
follows that the interface heights satisfy h(n + l,t) > 
h(n,t) + \ ——. Periodic boundary condition implies y(n± 
N) = y(n)±L, and correspondingly, h(n±N,t) — h(n,t). 

The dynamics of the interface involves the following 
moves in a small time interval dt: the move h(n, t) — > 
h(n,t + dt) = h(n,t) + 1 occurs with probability pdt 
while the move h(n,t) — > h(n,t + dt) = h(n,t) — 1 
takes place with probability qdt. The interface height 
remains the same with probability 1 — (p+ q)dt. The at- 
tempt to increase, respectively decrease, succeeds only if 
y(n+l,t) — y(n,t) > 1, respectively y(n,t) — y(n—l,t) > 
1. 



ASEP 


Lattice Interface Model 


Particle Label n 


Spatial coordinate n 


Displacement y(n, t) 


Height h{n, t) + 2 


Mean Velocity vp 


Mean growth rate ^ffi- 



TABLE I: Mapping ASEP to the interface. 



Following the above prescription, the ASEP with N 
particles and L sites maps onto a lattice model of the 
interface of length N . In order to get the continuum 
equation for the interface, we (i) coarse-grain the par- 
ticle labels so that the discrete tag label n becomes the 
continuous tag variable x, and (ii) divide x by the particle 
density p to make x into a spatial variable running be- 
tween and L. The equation of motion of the interface, 
to lowest order of nonlinearity, is given by 

dh(x,t) d 2 h dh A f dh\ 2 . . 



where vp = (p — q)(l — p),T = |,it = p(p — q),X — 
—2p(p— q). Here, r)(x,t) represents a Gaussian noise 
with (r)(x,i)) = 0, (rj{x, t)rj(x' , t')) = 2A5(x- x')5(t-t'), 

where A = | \ p)' boundary condition h(n ± 

N, t) = h(n, t) on the height variable for the discrete 
interface now reads h(x±L, t) = h(x, t) in the continuum. 
The derivation of Eq. [7] is relegated to the Appendix 
iBl The constant term vp and the drift term itlp in 
Eq. [7] can be eliminated by a boost and a Galilean shift, 
respectively, as explained below. In that case, Eq. [7] 
describes the time evolution equation of a KPZ interface 
[IH . Note that Eq. [7] is not an exact description of the 
time evolution for the ASEP density profile. It is rather 
a coarse-grained description of the ASEP. One expects 
that the scaling properties of the correlation functions for 
both the ASEP and its interface-equivalent are governed 
by the same KPZ fixed point, and hence, are described 
by the same critical exponents and scaling functions. 

For the symmetric exclusion process (SEP), the corre- 
sponding interface is an equilibrium one that does not 
move bodily. Rather, it fluctuates about a station- 
ary profile, and its time evolution is governed by the 
Edwards- Wilkinson (EW) equation [2(| , obtained by set- 
ting p — q, implying vp = 0, u = 0, A = 0, so that 



dh(x,t) d 2 h 



(8) 



On including all the higher order nonlinearities, Eq. [7| 
is replaced by the general form 



dh(x,t) _ T d 2 h 
dt dx 2 



m—0 



dx 



+ V(x,t), (9) 



where, e.g., Ao = vp, Ai = u, A 2 — A/2, etc. The different 
coefficients A m in Eq. [5] determine the scaling properties 
of the height fluctuations, as discussed below. Ao can be 
eliminated by redefining the height variable h — > h! = 
h + \ t. The first-order gradient term, Ai||, can be 
eliminated from Eq. [9]by Galilean shift x — > x' — x + \±t. 
The stationary state height fluctuations are measured 

by 



S(x,t) = ([h(x + x',t)-h(x',0)f) 



(h(x + x',t) - h(x',0)) 



(10) 



where the angular brackets denote averaging over the 
stochastic evolution of the interface, while the overbar 
represents averaging over the initial stationary ensemble 
of the interface configurations. 

If Ai is nonzero, the autocorrelation function grows as 
S(0,t) - t for large i 0,0,113 ■ This can be explained in 
terms of sliding density fluctuations (SDF) with respect 
to the particle motion, as discussed in Section IIVI If 
Ao = Ai = 0, the function S(x, t) assumes the scaling 
form 



(11) 
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in the asymptotic limit x,t — > oo, with x/t 1 ^ = con- 
stant. Here, (3 is the growth exponent while z is the dy- 
namical exponent. These exponents as well as the scaling 
function Y are the same for all systems belonging to the 
same universality class. The scaling function Y(s) has 
the property that Y(s) — ► constant as s — > 0, implying 
that the height autocorrelation S(0,t) ~ t 2/3 . Further, as 
s — > oo, the function Y(s) — > s 2a . Here, a is the critical 
exponent related to z and (3 through z = a/ (3 [2l|. It 
determines the roughness of the interface, S(x, 0) <~ x 2a . 

From a perturbative renormalization group (RG) anal- 
ysis of Eq. [9l it is known that there are two distinct fixed 
points in the parameter space. These two fixed points de- 
fine two universality classes. The first class is the one in 
which all A m = 0. In this case, Eq. [5] reduces to the EW 
equation, Eq. [8l The resulting linear growth problem 
was also studied by Hammersley in a different context 
[22| . The values of the critical exponents for this class 




= 1/4, z = 2 (EW in Id). (12) 

The Edwards- Wilkinson (EW) fixed point also controls 
the behavior in systems where only odd-order terms A m 
with m > 3 are present. The lowest-order nonlinearity 
is provided by the cubic term. Under a scaling trans- 
formation x — ► x' — bx, t — ► t' = b z t, h — ► h' = b a h , 
we see that the cubic nonlinear term is marginal around 
the EW fixed point. Calculations based on mode cou- 
pling theory [23[ and RG [2J, |25| show that it is actually 
marginally irrelevant, leading to multiplicative logarith- 
mic corrections to the power law: S(0,t) ~ i 1 / 2 (ln t) 1 / 4 , 
a behavior that has been verified by numerical simulation 

study mEim. 

Even-order coefficients lead to a breaking of h — ► — h 
symmetry of Eq. [HI causing a change in the universal- 
ity class. The second-order perturbation A2 is a relevant 
perturbation for the EW fixed point that drives the sys- 
tem away from the EW to a new KPZ fixed point. The 
KPZ fixed point is characterized by the exponent values 

m 

(3=1/3, z = 3/2 (KPZ in Id). (13) 

If Ai ^ 0, under the scaling transformation x — > x' = 
bx,t -> t' = b z t, ft, -> h! = b a h, both the EW and KPZ 
fixed points are unstable, and a flow towards a third fixed 
point (SDF) is generated (Fig. Q|, as we show in this 
paper. 



III. REPRESENTATION FOR AN INFINITE 
SYSTEM 

On the basis of the discussion in Section [Til wc con- 
clude that in the stationary state, the displacement of the 
n-th tagged particle, y(n,t) — y(n,0), has contributions 
from three distinct physical sources, namely, systematic 



SDF 

~ t 




(A m = V m) (Ai = 0) 

FIG. 4: Schematic representation of relative stabilities of var- 
ious fixed points, showing the associated behavior of the vari- 
ance <J 2 {L, i) of the displacement of the tagged particle in an 
infinite system. 



drift, sliding density fluctuations, and their dissipation. 
Correspondingly we write 

y(n, t) - y(n, 0) w v P t + t a G n (t) + t" X n(t), (14) 

where a = axpz = 1/2 is the roughness exponent while 
(3 = (3k pz = 1/3 is the growth exponent of the inter- 
face equivalent to the ASEP. The representation, Eq. [lj] 
is based on the following arguments. First, the coarse- 
grained density fluctuations for the ASEP are described 
by the KPZ interface whose fluctuation profile moves 
with velocity Vk ~ Vp. As a result, in time t, the particle 
senses the fluctuations over a stretch of length (vk — vp)t. 
From the scaling analysis of the interface summarized in 
Section|TTJ the typical fluctuation over this distance scales 
as [(vk — vp)t] a ~ t a . This fact is encoded in the second 
term in Eq |14l where G n (t) is a random variable which 
depends only on the initial configuration, drawn from the 
stationary ensemble, but is independent of the stochastic 
noise in the evolution of configurations. The fluctuations 
arising from the dissipation typically grow with time as 
f 3 and are represented by the last term in Eq. [141 Here 
Xn (t) is a random variable that depends only on the noisy 
history in the evolution of configurations, but is indepen- 
dent of the initial configuration. 

We show below how the form, Eq. [TU accounts for the 
observed behavior of correlation functions for an infinite 
system. 

• On averaging with respect to both the initial sta- 
tionary ensemble and stochastic evolution, we get 

(y(n, t) - y{n, 0)) w v P t + t a G^ + t { Xn ), (15) 

where, as previously, the angular brackets denote 
averaging with respect to the stochastic evolution, 
while the overbar is used to denote averaging with 
respect to the initial stationary ensemble. 
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It follows that for the variable defined in Eq. [TJ 
A„(i) « t a (G n -G^) + t { Xn - (Xn))- (16) 

With a = 1/2 and /? = 1/3 for the ASEP, the first 
term on the rhs dominates for large t. We recover 
the result for a 2 (L,t) = (A^(i)) for an infinite sys- 
tem (described in the Introduction), given by 

lim a 2 (L, t) 1 t(G n - G~ n f ~ t. (17) 

L — >oo 

• For fixed initial condition, averaging over noise his- 
tory, we get 

(y(n, t) ~ y(n, 0)) « « P t + TG„ + ^< Xn ). (18) 

Thus, for the displacement variable defined in Eq. 
02 we have 

dn{t)nt [ X n-{Xn)\. (19) 

The fluctuations in d n {t) are measured by s 2 (L,t) 
according to the definition, Eq. |U Thus, one has 

s 2 (t) = lim s 2 (L, t) ^ t^(( Xn - ( Xn )) 2 ) ~ (20) 

L — >oo 

With /3 = 1/3, we get the result s 2 {t) ~ i 2 / 3 Q 
for the growth of the variance of the displacement 
of a tagged particle in an infinite system, starting 
from an arbitrary fixed initial configuration drawn 
from the stationary ensemble, as discussed in the 
Introduction. 

IV. VARIANCES a' 2 (L,t) AND s 2 (L,t): 
PHYSICAL ARGUMENTS 

A. Infinite system 

The variance a 2 (t) of the displacement of a tagged par- 
ticle in an infinite system, averaged with respect to both 
the initial stationary ensemble and stochastic evolution, 
behaves very differently from the variance s 2 (t) where 
the averaging is with respect to only the stochastic evo- 
lution. Such differences for an infinite system have also 
been observed for shock fluctuations in the ASEP [27l ]. 
In this subsection, we provide physical interpretations of 
the observed linear growth in time for a 2 (t) and the t 2 / 3 
growth of s 2 (t). 

1. a\t) 

As the averaging in Eq. [2] is with respect to both an 
initial stationary ensemble and stochastic noise in the 
evolution of configurations, the average drift in time t is 
vpt, and fluctuations in the tagged particle displacement 
are defined with respect to this. In the rest frame of the 



density fluctuations, the tagged particle has an average 
velocity u = vp — In time t, it traverses a sequence 
of density fluctuations over a distance (up — vk)^, with 
each fluctuation adding a stochastic noise to the motion 
of the tagged particle. The noise is uncorrelated since the 
stationary state configurations have a product measure. 
The variance a 2 (L,t) is thus proportional to t for large 
t, by virtue of the central limit theorem. The coefficient 
of proportionality Do is known to be vp — (p — q)(l — p) 
0. This expression for D can be easily derived using 
the above picture of drift of the tagged particles relative 
to the density fluctuations [l2j]. 

2. s 2 (t) 

For s 2 (t), the fluctuations are measured by starting 
from a single fixed configuration, drawn from the ensem- 
ble of stationary states. Thus, in every measurement, 
the particle passes through the same sequence of density 
fluctuations. However, dissipation of the density profile 
is different from history to history. 

Define the deviation in the displacement of a tagged 
particle from vpt by 

D n (t) = y(n,t)-y{n,0)-v P t. (21) 

The dissipation in the density profile being different for 
different histories shows up in D n (t), as depicted by the 
gray trajectories in Fig. [5] Notice that the mean over 
histories, (D n (t)), is nonvanishing, as depicted by the 
black curve in Fig. [5] 




FIG. 5: The gray curves show the distance covered, D n (t), in 
time t about the mean vpt by the n-th tagged particle for 10 
Monte Carlo runs for a single fixed initial configuration, drawn 
from the stationary ensemble of the TASEP. The black curve 
shows the mean displacement, (D n (t)), obtained by averaging 
over 500 histories. Here, L = 10, 000, the particle density 
p = 0.25. 

Since the variance s 2 {t) of the fluctuation in Eq. 2] is 
measured with respect to vpt+ (D n (t)), it is able to sense 
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the dissipation of the density profile in time, without the 
effects of SDF. From Eq. Q31 with a = 1/2 and (3 = 1/3, 
we get D n {t) — t 1 / 2 G„(t) + t x ' 3 Xn(t)) on averaging over 
stochastic evolution, we get 



(D n (t)) =t 1 ' 2 G n (t)+t 1 ' 3 (xn(t)). 



(22) 



The first term in Eq. [22] dominates, and thus at large 
t, (D n (t)) is primarily determined by the pattern of spa- 
tial density fluctuations in the initial configuration. Ig- 
noring the dissipative component (the second term in 
Eq. Eg), we expect (D n (t)) to be given by D*(t) « 
(p — q) J* dt'{(l — p(x)) — (1 — p)}. Making a change 
of variable from t to x, and noting that the Jacobian of 
transformation is l/(vp — vk), we get 



(p-q) 

Vp — VK 



y(n,0) + (v P -v K )t 



y(n,0) 



dx{(l-p(x))-(l-p)}. 

(23) 

The upper limit of the integral incorporates the relative 
distance moved by the particle, (vp — vx)t, neglecting 
fluctuations of 0(t 1/2 ) coming from the second term in 
Eq. m 

In Fig. [HI we have compared (D n (t)) with D^(t) for 
a randomly chosen initial configuration. (D n (t)) was ob- 
tained from simulations, while D^(t) was obtained by 
integrating the initial density profile, following Eq. 
We see that there is good agreement between the signals 
for (D n (t)) and D'^(t). Fluctuations over a distance 5x 
dissipate in a time t ~ (5x) z with z — 3/2. Thus, with 
increasing t, spatial fluctuations in the initial profile are 
smeared out over larger distances, a fact which is borne 
out by Fig. [6] We have also measured the overlap func- 
tion for the signs of the two signals, (D n (t)) and £>*(£), 
for 10 different randomly chosen initial configurations at 
stationarity, and found a mean value 0.8, which indicates 
a fairly good degree of correlation between the two sig- 
nals. 



B. Finite system 



We now turn to the behavior of a 2 (L,t) and s 2 (L,t) 
on a finite lattice, and explain the differences in their 
behavior on the basis of SDF. 




10000 



FIG. 6: Comparison between (D n (t)) and D^(t) for a fixed 
initial stationary configuration of the TASEP. (D n (t)) was 
obtained from simulations, and involves averaging over 1000 
histories. D„(t) was obtained by integrating the initial den- 
sity profile according to Eq. [53] (see text). Here, L = 10, 000, 
while the density p = 0.25. 



The oscillations do not continue forever because of the 
cumulative effect of dissipation. An estimate for the time 
taken to completely dissipate an initial pattern of density 
fluctuations is T 2 ~ L z , where the dynamical exponent z 
is known to have the value 3/2 [28]. For times t ^S> T%, 
when the initial density profile has completely dissipated 
away, the fluctuations are entirely due to the motion of 
the center-of-mass: a 2 (L.t) ~ D(L)t with D(L) ~ -7=. 

The scaling of D(L) with the system size L is obtained 
by matching the behavior of a 2 (L,t) across T%. 

The lower envelope A 2 (L,t) of the oscillations in 
o 2 {L, t) is determined by the dissipation of the stationary 
state density profile, and can be studied by the method 
of sliding tags 0, @, QJ, H, SS HI H Hj . In order to 
monitor the dissipation of the moving density profile, we 
need to correlate, at different times, the location of two 
different particles. Since in the frame of the density fluc- 
tuations, particles move with velocity u = vp — vk, we 
need to examine the function 



K 2 {L,t) = {[y{n' 1 t)-y(nM 
where 



([y(n',t)-y(n,0)]) \ 
(24) 



2 (L,t) 



n — put. 



(25) 



In the rest frame of the density fluctuations, a tagged 
particle has velocity u = vp — vk and would take time 
T = L/u to return to its initial environment of density 
fluctuations. Thus the variance of the tagged particle 
displacement should increase from at t = 0, reach a 
maximum at around T/2, and come down to almost zero 
at t = T. The difference from zero at t — T is due to 
dissipation in the density profile arising from stochastic- 
ity in the dynamics. This scenario recurs at integral T's; 
the time period of oscillations is T = L/u = L/[(p— q)p] 



The tag shift put accounts for the relative motion of 
the particles and the density profile, and ensures that 
the time evolution of the same density patch is being 
recorded at every instant. Equationl25[ representing slid- 
ing of the tag, is tantamount to a Galilean shift that 
gets rid of the drift term, w§|, in Eq. \T\ In Fig 
we show how the method of sliding tags discussed above 
gives the lower envelope A 2 (L,t) of the oscillatory quan- 
tity a 2 (L,t). This envelope grows with time as t 2 ^ 3 un- 
til times - L 3/2 , beyond which A 2 (L,t) - D(L)t with 
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D(L) ~ -|=. Thus, beyond T 2 ~ L 3/2 , there is no dis- 
tinction between the temporal behavior of <j 2 (L,t) and 

a 2 cm). 



1200 



600 




20000 



FIG. 7: (Color online) MC simulation results for various 
tagged particle correlations. The topmost curve refers to the 
variance a 2 (L,t) of the displacement when the averaging is 
over both the initial stationary ensemble and stochastic evo- 
lution. The middle curve is the sliding tag correlation function 
A 2 (L,t) defined in the text; it coincides with a 2 (L,t) at the 
local minima of the latter. The lowermost curve shows the 
variance s 2 (L,t) when the averaging is only over stochastic 
evolution. Here, p — 0.25. The system size is 1024. The 
averaging is over 10 5 MC runs. 



2. 



'-(L,t) 



Since we always use the same initial condition in this 
case, the particle moves through the same sequence of 



density fluctuations in every measurement. Neverthe- 
less, the dissipation of the density profile is different for 
different histories, and s 2 (L,t) captures this. Typical 
fluctuation grows with time as i 1 / 3 , following Eq. Q31 
This leads to s 2 (L,t) ~ t 2 / 3 . This behavior continues 
until t T* ~ L 3 / 2 when the fluctuations are due to 
the diffusive motion of the center-of-mass. This gives 



! (L,t) - D(L)t with D(L) 



l 



One obtains the scal- 



ing of D(L) with the system size L by a simple match in 
the behavior of s 2 (L, t) across T* ~ L 3 ^ 2 . 



V. LINEAR INTERFACE MODEL: 
ANALYTICAL RESULTS 

The existence of size-dependent time scales in a 2 (L, t) 
and s 2 (L, t) for the ASEP, discussed in this paper, can be 
explained qualitatively by exactly solving the continuum 
equation, Eq. in the limit A — > 0. This is possible as, 
in this limit, Eq. [7] is linear, and hence, solvable. We 
refer to the corresponding interface as a linear interface. 
The relevant time evolution equation is 



dh(x, t) 
dt 



d 2 h dh , . 
Vp + T d^ +U ^ + ^- (26) 



In this section, we outline the exact computation of the 
two quantities <j 2 (L,t) and s 2 (L,t) for the linear inter- 
face. We will see that the drift term, in Eq. [2H that 
makes it different from the usual EW equation, plays a 
crucial role in determining the temporal behavior of these 
functions. 

Our aim is to compute a 2 (L,t) and s 2 (L,t) for the 
linear interface. For this, we look at the function 
CL(to,t + t), defined in Eq. [5J which, utilizing the map- 
ping outlined in Table[T]with n replaced by the continuous 
variable x. reads 



C L (to,t + t) = ([h{x,t + t) - h(x,t ) - ([h(x,t + t) - h(x,t Q )})} 2 ). (27) 

I 



As discussed in Section UH a 2 (L,t) can be obtained from 
the limiting behavior of the function Cl (to ,to + t) in the 
limit to — > oo. On the other hand, one recovers s 2 (L,t) 
from the limiting behavior of the function Cl^q, to + t) 
in the limit to — > 0. 

To compute Ci(io, t +t), we need to solve Eq. [26] for 
h(x,t). Before doing so, we note that the constant term 
on the right of Eq. [55] can be gotten rid of by going to a 
co-moving frame moving with velocity v p . This is equiv- 
alent to making the transformation h — * h + vpt. The re- 
sultant equation is the usual time-evolution equation for 
the EW interface, Eq. [5] with an additional drift term, 



u^, and can be solved by going to Fourier space. To this 
end, utilizing the boundary condition h(x±L, t) = h(x, t) 
on the height variable, we write h(x, t) in terms of its 
Fourier modes h(m,t). We have 

oo 

h(x,t)= h(m,t)e^ (x+ut) . (28) 

m— — oo 

Thus, h(m,t) = \ J Q L dx h(x, ^e-^^+^l . Equation 
|2"61 now reads 

dh(m,t) 47r 2 m 2 r , , _, , . . 

W =-r— =a-ftfrM)+jj(m,t), (29) 
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where 77(771, i) — j- j Q L dx rj(x,t)e t2 jf i i x + ut ) . Utilizing 
(r](x,t)r](x' ,t')) = 2A6(x - x')S{t - t'), we have, for the 
Fourier modes, 

2A 

(fj(m, t)f)(m', t')) = —A m ,_ m ,6(t - t'), (30) 

where A TOjn is the Kronecker delta. 

Equation [29] can be solved for h(m,t) to get 

h(m,t) = /i(to, 0)e~ l' 2 

_ r 4*rV a f f ,.i ~, r 4 ^ 2 " 2 *' 
+ e i- 2 at r)(m,t )e ^ 2 

Jo 

(31) 

The condition (rj(x, t)) = gives, for its Fourier modes, 



(77(771, t)) = V m. Using this in Eq. [31] we get 

(h(m,t)) = (h(m,0))e~ Ti ^ 2 ' t . 
From Eq. [3T1 with the help of Eq. [30] we get 



(h{m,t)h(m',t / )) = (h(m, 0)/i(m', 0))e ~ 



rl2_(m 2 t+m' 2 t') 



47r 2 r 



-r*^.|t-t'| _ e -r*^(t+f ) 



(32) 



A. <7 2 (L,t) 

Here, we need to take the limit to — ► oo in the function 
Ch{t(), to + t). Thus, we will compute 



a 2 (L,t)= lim C L (t ,t Q + t)= lim ([h{x,t Q +t)-h(x,t )-([h{x,t + t)-h(x 1 t )]}} 2 ). 



(33) 



r 



In the limit to — > oo, the quantity (/i(x, to + t) — M x ) to)) 
goes to zero. Thus, the function a 2 (L,t) reduces to 
a 2 (L,t) = lim io _ >ao ([ft(a;,to+*) - t )] 2 ). 

Utilizing Eq. [32] we get 



2A 

(h(x,t)h(x,t')) = (h 2 (0,0)) + —mxa(t,t') 



AL 
~47r 2 r 



E 



m— — oo.m^O 



L 

-r*&fi\t-t'\ 



x e L 



t(t-i') 



Utilizing this expression in Cx(to+t, to), and then taking 
the limit to — > oo, we finally get 



2 , r , 2A AL ^ 1 

L 7T^I 777^ 
m= 1 



1 — e j- 2 cos 



/ 2ixmu 



-t 



(34) 



(36) 

This is an exact formula for a 2 (L,t) within the linear 
model. Next, we consider the various limits. 

• t < L/u. Let 2 ^m = k. With Am = 1, we get 
Ak = 2^i. Now, ut < L implies that Afc is small, 
and hence, one can replace the sum over m in Eq. 
1361 by an integral over k to get 



In the last equation, the contribution from the initial con- 
dition, (h{m,Q)h(m' ,Q>))e- Ti ^ {mH+m ' H ' ] for m,mf + 
has been dropped, since, eventually when we set t = to 
and t' = to + 1 and let to — * °Oj this exponential term 
goes to 0. By combining the m and — m terms in the 
summation in Eq. 1341 we get 



a 2 (L,t) 



2Aut 



dk 



tiT Jo k 2 



cos(fc) 



(37) 



where, in obtaining the last equation, we have 
dropped the first term on the right of Eq. [36] in 
comparison with the second term. The integral on 
the right can be done exactly, see Appendix [C] Us- 
ing its value, we get 



2A 

{h{x,t)h{x,t')) = (h 2 (0,0)) + — min(M') 



A T OO - 

E^ 



27r 2 r ^ m 

rn—1 



/ 2itmu , ,. 
x cos I — - — (t — t) 



V L 



(35) 



a 2 (L,t) 



Au 



2—^Vter^ +t erf 



2%/T 



(38) 



Now we consider two cases: 



(a) t <C |£: In this limit, the error function is 
approximately zero, and we get 



a 2 (L,t) 



2A 



(39) 
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(b) « «t«£. Here, 



9 , T , Am 
<7 2 (L,f)« — i. 



(40) 



t - ^: Let f = ^ with n e I, Also, let 27r ^ = 
fe. Substituting in Eq. 1361 after replacing, for large 
L, the sum over m by an integral over k, we get 



L,t 



nL 



2At 



7T V Ml 



(i - e - fc ") 



( 41 ) 

The integral on the right can be done exactly, and 
its value is y/rm. This gives 



<r 2 (L,t 



nL, 
u 



2A 2A r 
L y/nT 



(42) 



For large L, keeping t fixed, when the first term on 
the rhs goes to zero, we have 



a 2 {L,t 



nL 



2A 



ttT 



V~t. 



(43) 



t^>L 2 : In this limit, the exponential term in the 
sum on the rhs of Eq. [35] drops out for all m to 

give 



9 2 A AL 

a 2 (L,t)n —t + — 

m— 1 



oo 1 



2A 



AL 

7r 2 r 



(44) 



In the limit of large t, the first term on the rhs 
dominates so that 



a 2 (L,t) 



2A 
T 



-t. 



(45) 



Note that the large time behavior of a 2 (L,t) is 
determined by the zero mode (m = 0) in the Fourier 
expansion of h(x,t). 

With the values of the exponents and z for EW 
given in Eq. I12[ we summarize the behavior of 
a 2 (L,t) for the linear interface in Table ITTl 

The scaling of the diffusion constant D(L) as the 
inverse of the system size in the regime t 3> T 2 ~ L z 
with z = zew — 2 can be obtained by matching the 
behavior of the a 2 (L,t) across the time T 2 . 

Let a 2 (L,t) ~ D(L)t with D(L) ~ L 7 for t > T 2 . 

On the lower side of T 2 , we have <r 2 (L, t) ~ t 2 @ . 

Then, to match the t and L dependence at T 2 , we 
must have 



{T 2 fP ~ D(L)T 2 . 



(46) 



This gives 2/3 = j/z + 1 whence, with z = 2 and 
(3=1/4 for the linear model, 7 = — 1. 



t < Ti ~ L: 
(a)*«£ • 
(b) ir «f«T 1 . 


a 2 (L,t) « 4?t. 


Ti < t < T 2 ~ L z . 

z = ZEW = 2. 
t = it with n e I. 

U 


„2/ r ,\ 2A .2/3 

/3 = /3ew = 1/4. 


t > T 2 . 


<T a (L,t)«^t. 
The diffusion constant D(L) = 2£. 



TABLE II: Behavior of a (L,t) in different time regimes for 
the linear interface. 



* < Ti ~ L: 
(a) t« f • 
(b) g <*<Ti. 


a 2 (L,t) ~ v 7 *- 
a 2 (L,t) ~ t. 


Ti < t < T 2 ~ L*. 

2 = z K pz = 3/2. 
t = — with tie/. 


a 2 (L,i) ~i 2/3 . 
(8 = /3xpz = 1/3. 


t > T 2 . 


~^t. 

The diffusion constant D(L) ~ 1/yL. 



TABLE III: Behavior of a 2 (L,t) in different time regimes for 
the KPZ interface. 



On the basis of the above results for the linear 
model, we expect the time-dependence of a 2 (L,t) 
for the KPZ class as in Table lTTll Here, z = zkpz = 
3/2 and = K pz = 1/3 so that 7 = -1/2. This 
value for 7 leads to the inverse square root scaling 
of the diffusion constant D(L) with system size L. 



B. s 2 {L,t) 

Here, we need to take to —> in the function Cx (to, £o + 
t). Thus, we will compute 

s 2 {L,t) = lim C L (t ,t +t) 
t — »o 

= ([h(x, t) - h(x, 0) - ([h(x, t) - h(x, 0)])] 2 ). 

(47) 

Now, from Eq. 1311 we have 



h(x, t) - h(x, 0) = [H m ,°) 



e ^ 2 + 



m— — 00 



S 1 -* / ^(m,t')e r ^-*'] xe^ 1 ^' 

OO 

/i(m,0)e l2 ^ a; . 



(48) 



Noting that every time we start from the same initial 
condition so that (h(m, 0)) = h(m, 0), we have 

h(x, t) - h(x, 0) - ([h(x, t) - h(x, 0)]) = 



12 



„- o Jo 



(49) 



Utilizing Eq. [49] and Eq. l30lm Eq. E3 we get, after a 
few steps, 



s 2 (L,t) 



oo 1 



AL 



47T 2 r * — ' TO" 1 



8 2 2 

e- r ^^*). (50) 



Separating out the m = mode and exploiting the m — > 
— m symmetry of the remaining terms, we get the final 
expression for s 2 (L, t) as 



9/ . 2A AL 

771=1 



OO ^ 



8 2 2 

e" r T^*). (51) 



Next, we will consider the various limits. 

• t < L 2 : Let fc = 2V2 [ 7rm v/t. Since Am = 1, we 
have Ak = ^^pLy/t is small for i < L 2 . This 
allows the sum over m in Eq. [51] to be replaced by 
an integral over k to give 



o, T , 2A i4 /2 r- Z" 00 dfc,., 
s 2 (L,t) = —t+-J-Vi r Tjl-e- 
Z, 7r V 1 1 k z 



-). (52) 



t < r* ~ ll 

2 = zew = 2. 


S 2 (L,i)«C7t 2/3 . 
/3 = Pew = 1/4. 

C = V^. 


t > r*. 

2 = ZEW = 2. 


The diffusion constant D(L) = 



TABLE IV: Behavior of s 2 (L,t) in different time regimes for 
the linear interface. 



Knowing the values of the exponents (3 and z for EW, 
given in Eq. 1121 we summarize the behavior of s 2 (L,t) 
for the linear interface in Table HVl The fact that D(L) 
scales as the inverse of the system size in the regime t 3> 
T* ~ L z with z = zew — 2, can be obtained by matching 
the behavior of s 2 (L, t) across the time T*. 

Thus, let s 2 (L, t) ~ D(L)t with D{L) ~ L 7 for <> T*. 

For i < T*, we have s 2 (L, i) - t 2f} . 

Then, to match the t and L dependence at T* , we must 
have 



(T*) 2/3 ~ D(L)T* 



(56) 



This gives 2/3 = 7/3+1. With /3 = 1/4 and z = 2 for 
the linear model, 7 = — 1. 

The above results lead us to expect the time- 
dependence of s 2 (L,t) for the KPZ class as in Table [V] 



Since t -C L 2 , the lower limit of the integral on the 
right can be taken to be 0. Then the integral can 
be done exactly, and its value is y/ir. For large L, 
keeping t fixed, the first term on the right of Eq. 
[55] goes to zero, and we get 



\L,t) w A\J^Vi. 



(53) 



• t ^> L 2 : In this limit, the exponential term in the 
sum on the rhs of Eq. [51] gets damped out for all 
to so that we finally have 



2fT . 2A AL ^ 1 
s 2 (L,i)«— i+ — £ — 

m — 1 



2A AL 
~L + 2tt 2 F V 6 



(54) 



For large i, the first term on the rhs dominates so 
that 



s 2 {L,t) w — t. 



(55) 



Here, as for the function <r 2 (L, t), we see that the 
large time behavior of s 2 (L, t) is determined by the 
zero mode (to = 0) in the Fourier expansion of 
h(x, t). 



t < T* ~ L z . 
z — zkpz = 3/2. 


s 2 (L,i) -i 2 ^. 
/3 = Pkpz = 1/3. 


t > r*. 


The diffusion constant D(L) ~ 1/y/L. 



TABLE V: Behavior of s 2 (L,t) in different time regimes for 
the KPZ interface. 



Here, (3 = (3k pz = 1/3 and z — zkpz — 3/2 so that 
2/3 = 7/2 + 1 gives 7 = —1/2. This explains the inverse 
square root scaling of the diffusion constant D(L) with 
system size L. 

Note that from Table [V] it follows that for an infinite 
system when the time-scale T* diverges, one recovers the 
result observed by van Beijeren that limL^oo s 2 (L,i) — 
s 2 (t)~t 2 / 3 i. 

To conclude this section, we reiterate that the ex- 
act solution for the linear interface qualitatively explains 
the occurrence of characteristic oscillations and different 
L— dependent regimes in the variance of the displacement 
on a finite lattice, as observed, for instance, in the the 
Monte Carlo simulations for the ASEP. We also note that 
in the large time limit it 3> L z with z = 3/2 for the 
ASEP and z = 2 for the linear interface), both <j 2 (L,t) 
and s 2 (L,t) grow linearly in time with the same con- 
stant of proportionality or the diffusion constant D(L) 
For the linear interface, D(L) ~ 
the ASEP. 



while D(L) ~ -7= for 
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VI. CORRESPONDENCE OF THE LINEAR 
MODEL TO DIFFERENT INTERACTING 
PARTICLE SYSTEMS 

The evolution equation, Eq. [26l in addition to being 
a linear approximation to the KPZ equation with a drift 
term, also arises in a number of microscopic interact- 
ing particle systems. Examples are the Katz-Lebowitz- 
Spohn (KLS) model at a specified value of the temper- 
ature, and the Asymmetric Random Average Process 
(ARAP). These correspondences are discussed below. 

A. Correspondence to the KLS model 

The aim of this subsection is to explain how the vari- 
ance of the displacement, computed in the linearized con- 
tinuum theory and defined in terms of the macroscopic 
parameters (T, A, speed u), compares to the same quan- 
tity evaluated by numerical simulations in the (discrete) 
KLS model, which uses one single microscopic parame- 
ter, namely, the temperature, for a special value of the 
latter. 



1. KLS model: a reminder 

The one-dimensional KLS model generalizes the ASEP, 
in that interactions between particles are added on top 
of the exclusion constraint [16| • The ASEP is the infinite 
temperature limit of the KLS model. 

The precise definition of the KLS model is as follows. 
Consider a chain of Ising spins s n (n = 1, . . . , A), evolv- 
ing under the Kawasaki dynamics with the heat-bath 
rule, and submitted to a drift. The energy of the chain 
reads 

E = - J)^s n s n+ x. (57) 

n 

In the heat-bath dynamics, for a pair of opposite spins 
(s„ + s n+1 = 0), the move (s„ -> -s n , s n+ i -> -s n+i ) is 
realized with probability 

W ( AE ) = 08) 

where AE — 2J(s„_iS n + s n +is n +2) = 0, ±4 J is the 
energy difference between the configurations after and 
before the move, and the numerator V is taken equal to 
p if the + spin is exchanged to the right (i.e., H — — > 

— h), and to q if it is exchanged to the left ( — I ► + — ). 

Associating a particle to a + spin, and a hole to a — spin, 
this particle hops to the right with probability p, and to 
the left with probability q. 

The moves corresponding to the three possible values 
of the difference AE are listed in the Table IVI1 with the 
corresponding acceptance probabilities. Evaporation cor- 
responds to the detachment of a 4- spin from a positive 



domain, or, equivalently, to the detachment of a parti- 
cle from a cluster of particles. Condensation, conversely, 
corresponds to the attachment of an isolated + spin to a 
positive domain, or, equivalently, to that of a particle to 
a cluster of particles, and the two diffusion mechanisms, 
to the motion of an isolated — spin in a positive domain 
(or hole in a cluster of particles) or to that of an isolated 
+ spin in a negative domain (or particle amongst empty 
sites). 

The heat-bath rule, Eq. [551 nas the non-trivial prop- 
erty that the steady state is independent of the asym- 
metry HIJ . It obeys detailed balance with respect to 
the energy, Eq. [57j at temperature T = l/(fcg/?) (&b: 
Boltzmann constant), i.e., W(AE) = W{~AE)e- 0AE . 



Type 


AE 


Acceptance Prob. 


Moves 


Condensation 


-47 


Pe 4/3J 7(e 4/3i + l) 


-+- + -*■ h + 

H 1 > + H 


Diffusion 





V/2 


++-+^+-++ 
h— <-» 1 


Evaporation 


+4J 


V/(e^ rj + 1) 


+ H ► + -+- 

h + ^- H h 



TABLE VI: Types of moves in the partially asymmetric 
Kawasaki dynamics, and corresponding acceptance probabil- 
ities with the heat-bath rule. The probability V is equal to 
p when the + spin exchanges to the right, to q when it ex- 
changes to the left. 



The KLS model can also be viewed as a migration 
process (or urn model). Generally speaking, a migra- 
tion process can be seen as a generalization of the Zero 
Range Process (ZRP) [H, |34|, where the particles hop 
on a lattice, but where the hopping rate depends both 
on the departure and arrival sites [35, 34]. Particles (or 
+ spins) are identified with the sites of a lattice, and the 
holes (or — spins) with particles located on these sites. At 
infinite temperature, this maps the ASEP onto the usual 
ZRP. At finite temperature, the rate at which a particle 
hops to the neighboring site is still given by Eq. [58] with 
AE = —AS J where 6 is the variation of the number of 
empty sites before and after moving the particle. This 
mapping is effectively used in numerical simulations. 

2. The variance of the displacement 

The displacement of a tagged particle in the KLS 
model obeys the same laws as for the ASEP, i.e., all 
phenomena presented in this paper for the ASEP exist 
likewise for the KLS model. More precisely, the tempo- 
ral behavior of the variance of the displacement of the 
tagged particle is the same as for the ASEP, up to pref- 
actors which depend continuously on the temperature. 
The question is now to make a link between the micro- 
scopic model and its continuum limit description, i.e., to 
compare the prediction of Section[V]for the temporal evo- 
lution of the variance, in the linear theory, to numerical 
simulations of the model. 
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For that purpose, we will (i) express the parameters 
of the continuum theory (r, A, and u) in terms of the 
temperature, which is the only parameter of the KLS 
model, (ii) explain the relevance of the linear theory for a 
special value of the temperature, for which the coefficient 
of the nonlinear term in the KPZ equation vanishes. This 
identification is possible in the present case because the 
stationary state of the KLS model is known. The method 
is given in [36| . It needs to be slightly generalized here, as 
now explained. For the sake of simplicity, we will restrict 
the computations to the case of zero magnetization M, 
or, equivalently, of density p — 1/2 (since M = 2p — 1). 

First, the following relationship between the diffusion 
constant T and the strength of the noise A, is always 
valid, independently of the asymmetry [36j ] . 



A = T X , 



(59) 



where \ is the susceptibility, explicitly known for the 
Ising chain, \ = e 2 ^ J . In order to find an expression of 
r and of the speed of the tagged particle, we write the 
equation for the temporal evolution of the magnetization 
(s„). This reads 



d(s n ) 
dt 



jsp 



jsp 

J n+li 



(60) 



where the spin current J* p through the link (n—l,n) has 
the following expression. 



40J 



2(eW + 1) 



((S„ - S n _l)(s„_ 2 - S„+l)) 



1 / e H ' J — 1 \ 

+ -((Sn-l) - (s n ) + , j -{(S n +l) - (s n -2) + (s„-2Sn-lSn) ~ \S n -XS n S n +l)) ) ■ 

-I - 2(e 4pj +1) 



(61) 



r 



The first line in Eq. [61] is the drift term. It is only 
present when the dynamics is asymmetric. It does not 
contribute to the computation of the diffusion constant 
r. Putting aside the drift term, the right-hand side of 
Eq. [HI] appears as a second order difference, consistently 
with the diffusive nature of the second term. (Note also 
that the first term is even in the spin variable, while the 
second one is odd.) In other words, the expression of T 
is the same both in the symmetric (EW) or asymmetric 
(KPZ) cases. It reads (see [36[ for the computation) 



1 



(eW + i)(eW + i) 



(62) 



It remains to evaluate the speed u, as well as the coef- 
ficient of nonlinearity A, appearing in the KPZ equation, 
Eq. [7] We use Eqns. IB11I and IB12I for this purpose. 
Now, without loss of generality, we consider completely 
asymmetric particle (i.e., + spin) motion with p — 1 and 
q = 0. For zero magnetization (or, equivalently, at den- 
sity p = 1/2), since u is given by u — vp 



\v P 



, dJ_i 
J dM\ 



M=0 



, while A is given by A 



M 
dp 

r> d 2 J 



p=l/2 



M=0 



need to evaluate the current of particles J as a function of 
the magnetization M. Note that J = J sp /2. We use Eq. 
I5T1 with (s n ) — M, (s n s n+ i) — (sxS2), by translation 
invariance, etc. Hence the current reads 



Note that this expression differs by a factor 2 with that 
given in [36[. This is due to the difference in the defini- 
tions of the scale of time chosen in the present work (see 
Table EU, and in (36|. Utilizing Eq. EH and the fact that 

the susceptibility x — e 2/3J , we get 



(63) 



( e -2/3J + 1 )( e 4/3J + 1 ) 



J =8 



1 



24/3,7 + 1 



{S1S3) 



2e^ J 
e 4 ? J + 1 



(sis 2 > 



(64) 



The correlators are known by the transfer matrix method 
which gives J as [36J 



2M 2 + e 4 ' 3j (l - M 2 ) - (1 + M 2 )(M 2 + e 4 ^ J (l - M 2 )) 1 / 2 
2( e 8/3J_i)(i_M 2 ) 

I 



(65) 



Utilizing the last equation, we finally get 



u = v P = 2J\ M =o = 



( e 2/3J + 1 )( e 4/3J + 1 ) ' 



(66) 



while 



A = 



1 - 3e 2 ^ J 



e 2/»J(eW + i)(e4W + i) 



(67) 



Hence for e 



-2/3 J 



3 (antifcrromagnetic chain), A van- 
ishes. For such a value of the "temperature" (actually of 
(3 J), though the dynamics is asymmetric, the continuum 
theory is linear. This yields 



r 27 9 a 9 

f = — , u = — , A = — . 

40 ' 40 40 



(68) 



Note that the unit of time is such that at infinite tem- 
perature 



r = 



A = 



l 



(69) 



These quantities differ by a factor 2 from their expression 
for the SEP (see Appendix [B| . As above, the origin of 
this difference lies in the definitions of the scale of time 
for these two models (see Table IVlf . 

The result of the comparison of the quantity a 2 (L,t) 
as obtained from simulation of the KLS model on a ring 
of size L for the temperature T satisfying e~ 2/3J = 3, 
and that from the exact solution of the linear interface 
as given in Eq. [35] (with the corresponding values of the 
constants F, A and u given in Eq. 155)1 is shown in Fig. 
[SJfor three different system sizes (L = 64, f 28, 256). The 
simulation of the KLS model was done in the equivalent 
migration process described above, with equal number 
of particles and sites, which corresponds to half filling 
(p = 1/2) in the original particle system (i.e. the KLS 
model). The variance of the displacement of a given 
tagged particle in the KLS model translates into the vari- 
ance of the number of particles which passed through a 
given bond in the migration process. 



B. Correspondence to the ARAP 

The Asymmetric Random Average Process involves 
hard core particles hopping on a real line as opposed 
to a lattice for the ASEP [Tj]. We consider a system of 
particles of average density p. We denote by Xi{t) the 
location of the z-th particle on the real line with peri- 
odic boundary conditions. The dynamics is stochastic 
and involves the following moves during an infinitesimal 
time dt: a randomly chosen particle jumps with proba- 
bility pdt to the right, with probability qdt to the left, 
and with probability f — (p + q)dt, it continues to occupy 
its original location. The amount by which the particle 
moves either to the right or to the left is a random frac- 
tion of the gap to the next particle to the right or to the 
left. Thus, for the i-th particle, the jump to the right is 
by the amount (xj+i — Xj), while to the left is by the 
amount by r~ (xj — Xj-i). Here, the random variables 
r i are independently drawn from the interval [0, 1], each 
being distributed according to the same pdf f(r), which 
is arbitrary. The time evolution equation of the positions 
Xi's is represented by the exact Langevin equation 




4000 



FIG. 8: (Color online) Comparison of the variance a 2 (L, t) of 
the displacement of the tagged particle in the KLS model ob- 
tained from simulations (dashed lines) and exact solution (full 
lines) of the corresponding continuum linear theory for the 
temperature T satisfying e~ 2J ^ knT = 3. Results are shown 
for three different system sizes (L = 64,128,256). Simula- 
tion of the KLS model is done with the equivalent migration 
process. 



where the random variables 7i(t) are given by 



7i(*) 



[Xi + i(t) ~ Xi(t)] with prob. pdt, 
\xi-i(t) — Xi(t)] with prob. qdt, 

with prob. 1 — (p + q)dt. 

The average velocity of a tagged particle vp is given by 



Vp 



d{x,{t)) 
dt 



<7i(*)>- 



(72) 



Xi(t + dt) =x i (t)+'Yi(t), 



(70) 



Using the expression for 7, in Eq. 1711 we can compute 
its average. This gives the average velocity of an ARAP 
particle as vp — iS—lltl^ where fx\ — f Q dr rf(r). Thus 
the average particle current J = pvp is independent of 
p. The derivation of the corresponding interface equation 
can be done along the lines outlined in Appendix [5) Not- 
ing that all derivatives of J with respect to p are zero, Eq. 
IB11I implies u — vp, while A = from Eq. IB121 Thus, 
we recover the equation for the linear interface, Eq. 1261 
Based on our exact solution for the linear interface, we 
expect the tagged-particle correlation in the ARAP in an 
infinite system as measured by <J 2 (t) = lim^oo a 2 (L,t) 
to vary as t (cf. Table HI)) . This has been confirmed by an 
exact solution of the tagged-particle correlation for the 
ARAP on an infinite system [37| . 
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VII. CENTER-OF-MASS MOTION 

In this section, we examine the motion of the center- 
of-mass for the ASEP, defined as 



1 N 



(73) 



n=l 



Define, for fixed initial configuration, drawn from the sta- 
tionary ensemble, 

d C M(t)=YcM(t)-YcM(0)-{YcM(t)-Yc M (0)), (74) 

where, as usual, angular brackets denote averaging over 
noise in the evolution of configurations. On the basis of 
the representation in Eq. I14| we have 

y(n, t) - y(n, 0) - (y(n, t) - y{n, 0)) = ^ [*„(*) - (*„(*)>]• 

This gives 



Hence, 



dcM{t) 



t 20 



N 



N 



£[*,(«) -(Xn(t))} 



(76) 



(4m W> =^E^(W')) - (Xn(t))(Xm(t))l 

. (77) 

The bracketed quantity on the right hand side, in the 
stationary state, becomes a function, f(n — m, t), of the 
difference. Putting this in Eq. [77] we get 



(d 2 CM (*)> 



~~N 



f(n - m,t). 



(78) 



The random variables Xn(t) are correlated only up to 
£(f) - t 1 / 2 . We write 

/•«(*) rl 
V f(n-m,t)=p dx f(x,t)~ pt llz \ dy g{y). 

n-m J ° J 

(79) 

In arriving at the second step in the above equation, we 
have made a transformation of the tag variable (n — m) 
to the spatial variable x. Plugging Eq. [75] in Eq. [75] we 
get 



(d 2 CM (t)} 



£2f3+l/z 



(80) 



t > L z , we expect {d 2 CM {t)) ~ D(L)t. Matching the be- 



havior of (d,Q M (t)) at the crossover time t* 
D(L) - L* 2 "- 1 ^. 

With z — zkpz = 3/2 and (3 = Pkpz 
ASEP, we have 



(d 2 CM (t)) 



ift«L 3 / 2 
* ift»L 3 / 2 . 



j L z , we get 
1/3 for the 



(81) 



The result (d 2 CM (t)} ~ ^- for * < L 3 / 2 for the ASEP 
has already been observed by van Beijeren et al. in [381 ] - 
For the linear interface in the EW class with j3 = 
Pew = 1/4, z = zew = 2, we get £>(L) ~ -p. Thus, 
for the linear interface, we expect, on the basis of our 
representation for the tagged particle displacement Eq. 




(d 2 CM (t)) 



L 



(82) 



for all times. Indeed, the above result, based on scaling 
arguments, is borne out by an exact computation of the 
quantity (d^ M (t)} for the linear interface model. Using 
the mapping of the ASEP to the interface in Table fl] for 
an interface of length L, we get 



d C M{t) = T / dx[h(x,t)-h(x,0)-([h(x,t)-h(x,0)})}. 
L Jo 

(83) 

We have h{m,t) = \ ^ dx h(x, t)e- l ^^ +ut \ This 

gives h(0, t) = jr J Q L dx h(x,t). Substituting in Eq. [83] 
we get 

dcM{t) = h(0,t)-h(0,0)-{[h(0,t)-h(0,0)]). (84) 
Using Eq. we finally get 



so that 



d C M{t)= / dt'f)(0,t') 
Jo 



(d 2 CM (t)) = f f dt'dt"(fj(0,t')fj(0,t")) 
Jo Jo 
2At 



(85) 



(86) 



where we have used Eq. [3D] in arriving at the last equa- 
tion. This exact result matches with the one obtained on 
the basis of scaling arguments in Eq. [521 



The above result is derived on the basis of the representa- 
tion for the tagged particle displacement, Eq. [14] based 
on a scaling argument. Using the mapping in Table[H Eq. 
1801 is expected to hold for an interface with nonequilib- 
rium dynamics with critical exponents [3 and z for times 
t <C L z . For large times t ^> L z , the fluctuations are 
expected to grow diffusively, as indicated by the linear 
behavior of the variances a (L,t) and s 2 (L,t). Thus, for 



VIII. CONCLUSION 

In this paper, we have revisited the problem of tagged 
particle correlation in the ASEP in one dimension, with 
particular emphasis on finite size effects. The stationary 
state variance of the displacement of a tagged particle 
involves averaging with respect to an initial stationary 
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ensemble and stochastic evolution. This quantity shows 
two distinct size-dependent time scales, T\ ~ L, and 
T 2 ~ L 3 / 2 . The behavior is linear for both the time 
regimes t -C T\ and t 3> T%, while for intermediate times 
T\ -C t -C T2, the variance <J 2 (L, t) shows pronounced os- 
cillations with a well-defined size-dependent time period. 
We understand the oscillations as arising from the slid- 
ing density fluctuations (SDF) relative to the drift of the 
tagged particle in the stationary state, the density fluctu- 
ations themselves being transported through the system 
by kinematic waves. 

Following van Beijeren we also study the variance 
s 2 (L, t) of the displacement of the tagged particle by aver- 
aging with respect to only stochastic evolution of a fixed 
initial configuration, drawn from the ensemble of station- 
ary states. We see that the quantity s 2 (L, t) has only one 
time scale T* ~ L 3 / 2 . The variance s 2 (L,t) grows lin- 
early in time for times t > T*, while for times t<T*, 
it grows as i 2 / 3 . 

Through a mapping to an interface, the time evolu- 
tion equation for the stationary state ASEP density pro- 
file maps onto that for a noncquilibrium interface in the 
Kardar-Parisi-Zhang (KPZ) universality class of inter- 
face dynamics. This equation is nonlinear and cannot 
be solved exactly. By dropping the nonlinear term, we 
obtain a linear model in the Edwards- Wilkinson (EW) 
universality class, which we solve exactly. This exact 
solution was helpful in understanding the occurrence of 
size-dependent time-scales Ti, T2 and T* , besides ex- 
plaining the scaling properties of two other interacting 
particle systems. 



Sp = p — (p), satisfies the equation of continuity 
d{5 P ) _ dj 



dt 



dx' 



(Al) 



where j is the particle current. Assuming that j depends 
on x through its dependence on a space-varying density 
p, the above equation can be recast into 



dp dp 

m = - c{p) dx-> 



(A2) 



where c{p) = j?-. To lowest order in Sp, Eq. IA2I gives 



dp 
dt 



TTT = —VK 



dp 

Ox 



(A3) 



with vk = c( (/?)). To leading order in Sp, we have 
Vk = |p) where J is the mean current in the steady 
state. Equation IA3I is the familiar wave equation in one 
dimension, and the solution reads Sp — g(x — vxt). Here, 
g is an arbitrary function of its argument. Thus, to lead- 
ing order in Sp, the solution to the equation of continuity 
gives kinematic waves of velocity vk that transport den- 
sity fluctuations through the system in the stationary 
state. The attribute 'kinematic' emphasizes the purely 
kinematic origin of these waves, in contrast to the dy- 
namic origin of say, elastic and acoustic waves, vk is 
identical to the collective velocity discussed in [Sd| ■ Note 
that this description neglects higher order terms in Sp, 
which is tantamount to neglecting dissipation in the den- 
sity profile while it is being bodily transported by the 
kinematic waves. 
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APPENDIX A: A PRIMER ON KINEMATIC 
WAVES 

The notion of kinematic waves goes back to the work 
of Lighthill and Whitham, who showed how their occur- 
rence follows from the continuity equation. We briefly 
recapitulate their argument [1, [39[ here. 

The steady state coarse-grained density fluctuation, 



APPENDIX B: DERIVATION OF THE 
INTERFACE EQUATION 

Here, we derive the interface equation, Eq. [3 following 
The local interparticle distance in the ASEP is given 

by 

1 = ^ (B1) 

p{x,t) dx' 

where y{x,t) is the position of the a;-th particle at time 
t. Now, 

y(x,t) = -+h{x,t), (B2) 

p 

where p is the mean density of particles and h(x, t) is 
the displacement of the ir-th particle from the position it 
would have had, had the particles been uniformly placed. 
Then, we have p(x,t)^ 1 = p^ 1 + implying 

KM) = A- ( B3 ) 

r Qx 
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Expanding in a power series in 



while 



p(x,t) 



,dh 



P 



Oh 



dp 2 



dx \dx 
The above equation can be rewritten as 



(B4) 



(B12) 



Using J(p) = (p — q)p(l — p) and vp = (p — q)(l — p) for 
the ASEP, we get 



dh 



p(x,t) = p+?p(x,t) where ip(x,t) — —p — — Vp' 

ox 



dh 

dx 



(B5) 

In the absence of any drift velocity, the equation of mo- 
tion of h(x, t) is diffusive. 



u = p(p- q), 
A = -2p(p-q). 



(B13) 
(B14) 



dh „,<3 2 /i . 



dt 



(B6) 



The noise term i](x,t) is Gaussian: (r)(x,t)) = 
0, (rj(x, t)r](x',t')) = 2AS(x - x')S(t -t'). 

In the presence of a drift velocity v(x, t), an additional 
term v(x,t) appears on the rhs of the above equation. 
The drift velocity depends on x and t only through the 
local density p(x, t) which is the only possibility in the 
coarse-grained lattice gas. Thus v(x,t) — v(p(x,t)). We 
now substitute p(x, t) = p+ip(x, t) and expand in a power 
series in tp(x, t). 



Utilizing the above expressions for u and A in Eq. IB10| 
we get the time evolution equation for the interface equiv- 
alent to the ASEP, Eq. El 

Expressions for the coefficients T and A in terms of 
microscopic parameters can be found by setting p = q = 
1/2, in which case the ASEP reduces to the SEP. The 
coefficient T can then be calculated explicitly [Z^, EH, 
with the result 



(B15) 



v(x,t) = v(p) 
Thus, we get 



dv 

d~ P 



ip(x,t) 



dp 2 



iP 2 (x 1 t) + 



(B7) 



Further, Eq. IB10I reduces to the EW equation, Eq. [5J 
and cr 2 (t) can be found from Section|V]by considering the 
limit u — ► 0. The result is <r 2 (t) ~ -^=y/t. Comparing 

with the exact result a 2 (t) w (^) s/i for the SEP 

H, one gets -fe = ^ 
get 



p 



dh 
~di 



V{p)+T dx^- P 



dv 


dh 


1 


d 2 v~ 




(dh\ 


dp 


p dx 


" 2 


dp 2 . 


P A 

p 


\dx) 



Using Eq. IB15) we finally 
(B16) 



dv 




fdh\ 


dp. 


P 3 

p 


Kdx) 



+ ... + r](x,t) 



, . w<9 2 /i .dh , , , 



~2 



dh 

dx 



(B8) 

+ T](x,t). 

(B9) 



APPENDIX C: EVALUATION OF THE 
INTEGRAL - e~ cfe2 cos(fc)] 



Let 



Here, v(p) is the mean drift velocity of the particle 
and hence, equals vp in our notation. Also, v! — 

v (p) ~ 7T~ where J(p) = pv(p) is the mean cur- 

rent. In terms of the kinematic wave speed Vk, we 
have v! = p(vp — vk)- The nonlinearity coefficient 

Note that x in the above equation stands 



r°° dk 

1(c) = -J[l-e-«* cos(fc)]. (CI) 
Jk=o K 



Thus, 



dJ 
dc 



dk e cos(k) = -\ ~e 
fc =o 2V c 



(C2) 



A' = p 3 

for the tag variable in the continuum. Now, dividing x by 
the particle density p to make it into a spatial variable, 
we finally get from Eq. IB9[ to lowest order of nonlinear- 
ity, 



Also, 1(0) = f . Hence, 



dh _ T d 2 h 
dt Vp dx 2 

where 



dh A / dh s . 

"s + ns,i + *' t '- (B10) 



'M-HWi'- 1 ' 4 '- (c3) 

Doing the integral on the rhs by parts, we finally get 



vp - v K , 



(Bll) 



1(c) = - + V^ z~ 1/ic ~ V / d V e~ V y- 1/2 - (C4) 

Z Z Jl/Ac 

Using the usual definition of the complementary error 
function, erfc(z) = ^= dt e~* = 1 — erf(z), (where 
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erf(z) is the usual error function), we can rewrite the 
above expression as 



7TC e 



-l /Ac 



7T 



- erf ( J, 

2 



— + V 7TC e 



"^-5"* (55s) (C5 » 



[1] T. M. Liggett, Interacting Particle Systems (Springer- [20 

Verlag, New York, 1985). 
[2] A. De Masi and P. A. Ferrari, J. Stat. Phys. 38, 603 [21 

(1985); R. Kutner and H. van Beijeren, J. Stat. Phys. 

39, 317 (1985). 

[3] B. Derrida, M. R. Evans, and D. Mukamel, J. Phys. A: [22 
Math. Gen. 26, 4911 (1993); B. Derrida and K. Mallick, 
J. Phys. A: Math. Gen. 30. 1031 (1997). 

[4] H. van Beijeren, J. Stat. Phys., 63, 47 (1991). 

[5] S. N. Majumdar and M. Barma, Phys. Rev. B 44, 5306 [23 
(1991). 

[6] S. N. Majumdar and M. Barma, Physica A 177, 366 

(1991). [24 

[7] In the absence of the external drive, it is known that 

a 2 {L,t) ~ t 1 ' 2 [T. E. Harris, J. Appl. Probab. 2, 323 [25 
(1965); R. Arratia, Ann. Probab. 11, 362 (1983); see 
also Ref. 0] ]. In the ASEP, if the external drive is small, [26 
there is a size-independent short-time precursor to the 
initial linear regime where a 2 (L,t) ~ t 1 ^ 2 . [27 

[8] M. J. Lighthill and G. B. Whitham, Proc. R. Soc. London 

A 229, 281 (1955). [28 

[9] M. J. Lighthill and G. B. Whitham, Proc. R. Soc. London 

A 229, 317 (1955). [29 
[10] D. Chowdhury, L. Santen, and A. Schadschneider, Phys. 

Rep. 329, 199 (2000) Section 4.1. [30 
[11] J. Lee, Phys. Rev. E 49, 281 (1994). 

[12] M. Barma, J. Phys. A: Math. Gen. 25, L693 (1992). [31 

[13] M. Barma, in Non Linear Phenomena in Materials Sci- 
ence III - Instabilities and Patterning, edited by G. Anan- [32 
thakrishna, L. P. Kubin, and G. Martin (Scitec Publica- [33 
tions, Untermiili, 1995). 

[14] M. Barma and R. Ramaswamy, in Non-Linearity and [34 
Breakdown in Soft Condensed Matter, edited by K. K. 
Bardhan, B. K. Chakrabarti, and A. Hansen (Springer- [35 
Verlag, Berlin, 1994). 

[15] M. Kardar, G. Parisi, and Y. C. Zhang, Phys. Rev. Lett. [36 
56, 889 (1986). 

[16] S. Katz, J. L. Lebowitz, and H. Spohn, Phys. Rev. B 28, [37 
1655 (1983); (b) J. Stat. Phys. 34, 497 (1984). 

[17] J. Krug and J. Garcia, J. Stat. Phys. 99, 31 (2000); [38 
R. Rajesh and S. N. Majumdar, J. Stat. Phys. 99, 943 
(2000). [39 

[18] D. Kim, Phys. Rev. E 52, 3512 (1995); O. Golinelli and 

K. Mallick, J. Phys. A: Math. Gen. 38, 1419 (2005). [40 

[19] Note that there are other mappings of ASEP to an inter- 
face that have their own merits, e.g., Ref. [1J; M. Plischke, 
Z. Racz, and D. Liu, Phys. Rev. B 35, 3485 (1987). Here, [41 
one identifies the presence of a particle with a downward 
slope of the interface while the absence of a particle is [42 
identified with an upward slope. In this way, the ASEP [43 
density profile maps onto a nonequilibrium interface, de- 
scribed by the KPZ equation. 



S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. Lon- 
don, Ser A 381, 17 (1982). 

A. -L. Barabasi and H. E. Stanley, Fractal Concepts 
in Surface Growth (Cambridge University Press, Cam- 
bridge, 1995). 

J. M. Hammersley, in Proceedings of the fifth Berkeley 
Symposium on Mathematical Statistics and Probability, 
edited by L. M. Le Cam and J. Neyman (University of 
California Press, Berkeley, 1967). 

B. Derrida, J. L. Lebowitz, E. R. Speer, and H. Spohn, 
Phys. Rev. Lett. 67, 165 (1991); (b) J. Phys. A. 24, 4805 

(1991) . 

M. Paczuski, M. Barma, S. N. Majumdar, and T. Hwa, 
Phys. Rev. Lett. 69, 2735 (1992). 

P. M. Binder, M. Paczuski and M. Barma, Phys. Rev. E 
49, 1174 (1994). 

B. Subramanian, G. T. Barkema, J. L. Lebowitz and E. 
Speer, J. Phys. A: Math. Gen. 29, 7475 (1996). 

F. J. Alexander, S. A. Janowsky, J. L. Lebowitz, and H. 
van Beijeren, Phys. Rev. E 47, 403 (1993). 

D. Dhar, Phase Transitions 9, 51 (1987); L. -H. Gwa and 
H. Spohn, Phys. Rev. A 46, 844 (1992). 

G. Tripathy and M. Barma, Phys. Rev. Lett. 78, 3039 
(1997). 

D. Das, A. Basu, M. Barma and S. Ramaswamy, Phys. 
Rev. E 64, 021402 (2001). 

A. A. Ferreira and F. C. Alcaraz, Phys. Rev. E 65, 052102 
(2002). 

J. M. Luck, C. Godreche, J. Stat. Mech. P08009 (2006). 
M. R. Evans and T. Hanney, J. Phys. A: Math. Gen. 38 
R195 (2005). 

C. Godreche, |cond-m at/0603249 also, Lecture Notes in 
Physics 716, (Springer- Verlag, 2007). 

F. Kelly, Reversibility and Stochastic Networks (Wiley, 
Chichester, 1979). 

C. Godreche and J. M. Luck, J. Phys. A: Math. Gen. 36, 
9973 (2003). 

R. Rajesh and S. N. Majumdar, Phys. Rev. E 64, 036103 
(2001). 

H. van Beijeren, R. Kutner, and H. Spohn, Phys. Rev. 
Lett. 54, 2026 (1985). 

G. B. Whitham, Linear and Nonlinear Waves (John Wi- 
ley & Sons, New York, 1974). 

G. M. Schiitz in Phase Transitions and Critical Phenom- 
ena 19, edited by C. Domb and J. L. Lebowitz (Aca- 
demic, London, 2001) Section 4.2.2. 

S. N. Majumdar, Ph. D. Thesis (Bombay University) 

(1992) (unpublished). 

R. Kutner, Phys. Lett. 81 A, 239 (1981). 

R. B. Stinchcombe, M. D. Grynberg, and M. Barma, 

Phys. Rev. E 47, 4018 (1993). 



